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^ 1 ABSTRACT 

q ■ We extend the formalism for the calculation of the relativistic corrections to the 

, Sunyaev-Zel'dovich effect for clusters of galaxies and include the multiple scatter- 

ing effects in the isotropic approximation. We present the results of the calculations 



■ by the Fokker- Planck expansion method as well as by the direct numerical integration 

of the collision term of the Boltzmann equation. The multiple scattering contribution 
is found to be very small compared with the single scattering contribution. For high- 
temperature galaxy clusters of ksT e w 15keV, the ratio of the both contributions is 
—0.2% in the Wien region. In the Rayleigh- Jeans region the ratio is —0.03%. There- 
fore the multiple scattering contribution is safely neglected for the observed galaxy 
clusters. 
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1 INTRODUCTION 

Compton scattering of the cosmic microwave background (CMB) radiation by hot intracluster gas — the Sunyaev-Zel'dovich 
effect (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972, 1980a, 1980b, 1981) — provides a useful method to measure the 
Hubble constant Ho (Gunn 1978; Silk & White 1978; Birkinshaw 1979; Cavaliere, Danese & De Zotti 1979; Birkinshaw, Hughes 
& Arnaud 1991; Birkinshaw & Hughes 1994; Myers et al. 1995; Herbig et al. 1995; Jones 1995; Markevitch et al. 1996; Holzapfel 
et al. 1997; Hughes & Birkinshaw 1998; Furuzawa et al. 1998; Komatsu et al. 1999; Reese et al. 2000). The original Sunyaev- 
Zel'dovich formula has been derived from a kinetic equation for the photon distribution function taking into account the 
Compton scattering by electrons: the Kompaneets equation (Kompaneets 1957; Weymann 1965). The original Kompaneets 
equation has been derived with a nonrelativistic approximation for the electron. However, recent X-ray observations have 
revealed the existence of many high-temperature galaxy clusters (David et al. 1993; Arnaud et al. 1994; Markevitch et al. 
1994; Mushotzky & Scharf 1997; Markevitch 1998). In particular, Tucker et al. (1998) reported the discovery of a galaxy 
cluster with the electron temperature fcsT e = 17.4 ± 2.5 keV. Rephaeli and his collaborator (Rephaeli 1995; Rephaeli & 
Yankovitch 1997) have emphasized the need to take into account the relativistic corrections to the Sunyaev-Zel'dovich effect 
for clusters of galaxies. 

In recent years remarkable progress has been achieved in the theoretical studies of the relativistic corrections to the 
Sunyaev-Zel'dovich effects for clusters of galaxies. Stebbins (1997) generalized the Kompaneets equation. Itoh, Kohyama & 
Nozawa (1998) have adopted a relativistically covariant formalism to describe the Compton scattering process (Berestetskii, 
Lifshitz & Pitaevskii 1982; Buchler & Yueh 1976), thereby obtaining higher-order relativistic corrections to the thermal 
Sunyaev-Zel'dovich effect in the form of the Fokker-Planck expansion. In their derivation, the scheme to conserve the photon 
number at every stage of the expansion which has been proposed by Challinor & Lasenby (1998) played an essential role. 
The results of Challinor & Lasenby (1998) are in agreement with those of Itoh et al. (1998). The latter results include 
higher-order expansions. Itoh et al. (1998) have also calculated the collision integral of the Boltzmann equation numerically 
and have compared the results with those obtained by the Fokker-Planck expansion method. They have confirmed that the 
Fokker-Planck expansion method gives an excellent result for fcsT e < 15keV, where T e is the electron temperature. For 
fcflTe > 15keV, however, the Fokker-Planck expansion results show nonnegligible deviations from the results obtained by the 
numerical integration of the collision term of the Boltzmann equation. 

In our previous papers devoted to the study of the relativisitc corrections to the Sunyaev-Zel'dovich effect for clusters of 
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galaxies (Itoh, Kohyama, & Nozawa 1998; Nozawa, Itoh, & Kohyama 1998; Itoh, Nozawa, & Kohyama 2000; Nozawa et al. 
2000), we have so far restricted ourselves to the case of single Compton scattering. This is justified because the optical depth 
for the Compton scattering of the CMB photon inside the galaxy clusters is generally about 10~ 2 or smaller (Birkinshaw 
1999). Nevertheless, it would be desirable to evaluate the effects of the multiple Compton scattering of the CMB photon 
inside the galaxy clusters accurately, as we have already developed the method to calculate the relativistic corrections to 
the Sunyaev-Zel'dovich effect for the galaxy clusters with high accuracy. The multiple scattering effects have been already 
considered by many authors (Wright 1979; Fabbri 1981; Loeb, McKee & Lahav 1991; Sazonov & Sunyaev 1998; Molnar & 
Birkinshaw 1999; see Birkinshaw 1999 for other references). Molnar & Birkinshaw (1999), in particular, have carried out a 
detailed Monte Carlo calculation including multiple scattering effects. However, most of the calculations to date other than 
Monte Carlo type have assumed isotropy of the radiation field after the first Compton scattering. In this paper we wish 
to evaluate the multiple scattering effects in the same theoretical framework of our previous papers. As a matter of fact, 
the method of calculating the multiple scattering contributions to the Sunyaev-Zel'dovich effect adopted by Fabbri (1981) 
and also by Sazonov & Sunyaev (1998) is in the same line as the present paper. The lowest-order term has been already 
obtained by them. We will calculate the relativistic corrections (higher-order terms) in the present paper. We will also carry 
out direct numerical integration of the collision term of the Boltzmann equation and compare the results with those obtained 
by the Fokker-Plank expansion method. The present paper is complementary to the work of Molnar & Birkinshaw (1999) 
who presented Monte Carlo results without the isotropic approximation. In the central region of a spherical galaxy cluster, 
the assumption of the isotropy of the incident radiation field is valid. Therefore, rigorously speaking, the present result is 
valid for such conditions of restricted geometry. Nevertherless, the analytical as well as numerical results with the isotropic 
approximation in the present paper will be useful when one compares with the full numerical calculation which does not 
assume the isotropy of the incident radiation field. 

The present paper is organized as follows. In § 2 we give the method of the calculation and the results. In § 3 we give 
discussion of the results and concluding remarks. 



2 MUTIPLE SCATTERING CONTRIBUTION IN THE ISOTROPIC APPROXIMATION 

In the present paper, we would like to derive the analytic as well as numerical expressions for the multiple scattering con- 
tribution to the Sunyaev-Zeldovich effect for the central region of a spherical galaxy cluster in the isotropic approximation. 
As a reference system, we choose the system that is fixed to the center of mass of the galaxy cluster. The galaxy cluster is 
assumed to be fixed to the cosmic microwave background (CMB). Following Itoh et al. (1998), we start with the Fokker-Planck 
expansion for the time evolution equation of the CMB photon distribution function n(o>): 
dn(u>) „ f dn ,„ .1 r 

et = 2 [d- x +n{1 + n) \ h 
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In equation (4), W is the transition probability of the Compton scattering, f(E) is the relativistic Maxwellian distribution 
function for electrons with temperature T e . We have integrated equation (4) analytically with power series expansions of the 
integrand. The expansion parameter is 
„ ksT e 



(5) 



The explicit forms for Ik are given in Itoh et al. (1998). 

We first assume the initial photon distribution of the CMB radiation to be Planckian with temperature To: 
1 



n(X) = n (X) 
where 
X 



z x -1 ' 



fcsT 



(6) 



(J) 



Assuming To/T e -C 1, one obtains the following expression for the fractional distortion of the photon spectrum derived by 
Itoh et al. (1998): 



An(X) = yJ^Ce? r q2 f qA q5 q6 , 



II 



G T / d£lJV e (€l), 



(8) 



(9) 



where ot is the Thomson scattering cross section, N e is the electron number density, and the integral is over the photon path 
length in the cluster. The explicit forms for Yo, Yi, Yi, Yz and Y4 are given in Itoh et al. (1998). In the present paper we have 
also included higher order terms of O^Ys and &\Y§. The explicit forms are given as follows: 
v 45 7515 ~ 28917 ^ 2 795429 ^ 3 2319993 ^ 4 12667283 ^ 5 

Ys - ~Y~^2T X ~^ X + ^^ X TT~ X + 112 x 

806524 - 6 21310 -7 46679 ^ , 10853 ^ 58 , 29 - u 
2i~ X + ^~ X ~-63~ X + ^52T X ~ 45 X + 1890* 

52 / 28917 795429 ~ 25519923 ^ 2 164674679 ^ 3 7661978 ^ 4 
t ij I ; — H ; — — A H z~rz -A — — Jv 



2N 



112 
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+ 4262oox 5 - 11529713 + 2 -^r - + ^r- 
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V ~ + 315 * ~ 225 X + 6300 * J 
( 1616456 34394053 ~\ 
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X = Xcoth(^) , (12) 

§ s — 7TV- (13) 

Equation (8) is the single scattering contribution, i.e. the first-order term in y. If the cluster of galaxies is optically thin, 
i.e. j < 1, the single scattering approximation is a good approximation. In fact, the approximation is valid for most of 
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the clusters. However, it is extremely important to calculate the next-order contribution in order to obtain more accurate 
theoretical prediction for the future observation of the Sunyaev-Zeldovich effect for clusters of galaxies. 

We now calculate the multiple scattering contribution. Since y <C 1 is realized for most of clusters of galaxies, the second- 
order contribution is considered to be sufficient. We now assume that the initial photon distribution has an isotropic first-order 
perturbation. Namely, 



n(X) = m(X) = n {X)+An(X), 

' An(X) 1 

^ iX) < 1 + ^xf\ 



(14) 



where the second term in equation (14) is given by equation (8). Here we have followed Fabbri (1981) and Sazonov & 
Sunyaev (1998) in assuming that the radiation field after the first Compton scattering is isotropic. To go beyond the isotropic 
approximation would be a very involved calculation depending on the exact geometry of the galaxy cluster. In such a calculation 
simple analytic expressions resulting from the Fokker-Planck expansions would not be possible. In this paper, therefore, we 
content ourselves with the isotropic approximation and carry out higher-order Fokker-Planck expansions. Inserting equation 
(14) into RHS of equation (1), and performing the standard calculation, we obtain the following expression for the fractional 
distortion of the photon distribution function including the second-order contribution: 
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=,10 / 4092665438636 22582640792464 - 8264177610763 ^ 2 180539814396049 ~. 3 
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451257055744 - 4 271877622784 
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where X and S are defined by equations (12) and (13), respectively. In equation (15), the first term corresponds to the first- 
odrder contribution and the second term corresponds to the second-order contribution of the multiple scattering. In deriving 
equation (15), we have used the following identity relation for y: 

o T dliN e {h)a T 1 de 2 N e (£ 2 ) = i (a T J dtiN e {£i)j 

= \v 2 - (23) 

We have also neglected terms higher than 0(0® ) in the y 2 contributions in equation (15). It is important to note that equation 
(15) satisfies the photon number conservation. 

We note that the lowest-order term Zo has been already obtained by Fabbri (1981) and also by Sazonov & Sunyaev 
(1998). The present calculation has produced the relativistic correction terms Zi, Z2, Z3, Z4, and Ze- 

We will also calculate the multiple scattering effect to the thermal Sunyaev-Zeldovich effect by numerically integrating 
the collision term of the Boltzmann equation. Here we assume that the radiation field after the first Compton scattering is 
isotropic. The original Boltzmann equation for the photon distribution function n(uj) is given by (Itoh et al. 1998) 

^p = -2 / ^d 3 p'd 3 k'Wf{E) {[l + n (,>H-[l + n H]n(,') e 41 } . (24) 

By inserting equation (14) in equation (24) we obtain 

dU{Uj) ~ -2/ ^dyd\'Wf(E){[l + n (Lo')]no(u)-[l + no(Lu)]no(Lo')e Ax } 



dt J (2tt 

-2 



^d 3 p'd 3 k' W f{E) L {u')^l [noH (1 - e *») - e A * 



+n (to)^± [1 + no(a /) (! _ e A*)] ) . (25) 
n (w) J 

The first term corresponds to the single scattering contribution which has been calculated by Itoh et al. (1998). The second 
term corresponds to the double scattering contribution. For An(u/)/no(w') we use the numerical results obtained by Itoh et 
al. (1998) and calculate the double scattering contribution by numerical integration of the collision term of the Boltzmann 
equation. 

With equation (15), we define the distortion of the spectral intensity as follows: 

Al = ^T^-A/ 1+ A/ 2 . (26) 
e x - 1 n (X) 
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The first term AIi contains a factor y, whereas the second term AI2 contains a factor y 2 . In Figure 1 we show A/2/y 2 as a 
function of X for the case fcsT e = 5keV. For fcsT e < 5keV it is found that results of the Fokker-Planck expansion approximation 
perfectly agree with that obtained by the numerical integration of the collision term of the Boltzmann equation for the entire 
region of X < 20. For ksT e < 5keV the convergence of the series with respect to the relativistic temperature parameter 9 e 
is relatively fast. Therefore it is sufficient to include up to B\Z/± terms in equation (15). In Figure 2 we show A/2/2/ 2 f° r the 
case fcflTj = lOkeV. For this temperature region the convergence of the series expansion is slow for large values of X. The 
Fokker-Planck expansion aprroximation is valid for X < 10. For higher temperature region, the convergence is even worse. 
For fcflT e = 15keV the expansion approximation is valid for the region of X < 4. 

In order to estimate the relative importance of the multiple scattering contribution, we now define the following ratio: 

r ■ 

In Figure 3 we show V for X — 5 as a function of the electron temperature T e . It is clear that F increases with a negative 
sign as the temperature of the cluster of galaxies increases. The numerical result (solid curve) shows that the maximum 
contribution is V w —0.2 at fesT e = 15keV. However, the multiple scattering contribution is small because of a further factor 
y. Namely, for the cluster of galaxies of fcsT e = 15keV, we have 

^ = yr «-0.2 V w-0.2%, (28) 

where we used a typical value y ~ 0.01 of the galaxy clusters. Therefore the maximum effect of the multiple scattering 
contribution is —0.2% of the single scattering contribution for the observed high-temperature galaxy clusters. 

In the Rayleigh-Jeans limit where X — * 0, equation (15) is further simplified: 
An(X) _ _ 17 0e+ 123^ _ 1989^ + 14403,4 _ 20157^ 423951^ 



n (X) a \ 10 40 e 280 e 640 e 224 e 1024 

+ 2y 2 ei fi - l le e + ™el - ^el + - ^^el + 25927827 el) . (29) 

y \ 5 25 1400 175 4480 19600 J v ' 

With equation (29) we have the multiple scattering contribution for fcgT e =15keV and y — 0.01 as follows: 

-Jl « -y0 e w -0.03%. (30) 

In the Rayleigh-Jeans region the multiple scattering contribution is safely neglected. 



3 DISCUSSIONS AND CONCLUDING REMARKS 

We have calculated the relativistic corrections to the multiple scattering contribution to the Sunyaev-Zel'dovich effect in the 
isotropic approximation by extending the formalism developed in our previous papers as well as by Fabbri (1981) and Sazonov 
& Sunyaev (1998). We have also calculated the multiple scattering effect to the Sunyaev-Zeldovich effect by numerically 
integrating the collision term of the Boltzmann equation. Our approach is complementary to the Monte Carlo calculation of 
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Molnar & Birkinshaw (1999). We have estimated the accuracy of the Fokker-Planck expansion approximation by comparing 
with the result obtained by the numerical integration of the collision term of the Boltzmann equation. We have found that 
Fokker-Planck expansion approximation is valid for the enire region of X < 20 for kBT e < 5keV. However, for fcsT e < lOkeV, 
valid region is limited to X < 10. For higher temperature region, the convergence is even worse. For ksT c = 15keV the 
expansion approximation is valid for the region of X < 4. 

From the results presented in the previous section it is clear that the multiple scattering contribution A/2 is very small 
compared with the single scattering contribution Aii. For high-temperature galaxy clusters of ksT e « 15keV, we obtain the 
ratio A/2/A/1 ~ —0.2% at X = 5. In the Rayleigh-Jeans region we have AI2/AI1 ~ —0.03%. Therefore it is concluded that 
the multiple scattering contribution to the thermal Sunyaev-Zel'dovich effect for galaxy clusters can be safely neglected. The 
reader is therefore referred to our previous four papers which deal with the single scattering contribution in detail. 
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FIGURE CAPTIONS 

• Figure 1. The spectral intensity distortion AI2/1/ 2 as a function of X for the case ksT e = 5keV. The dotted curve shows 
the contribution up to d\Z^. The dashed curve shows the contribution up to Q\Z$. The dash-dotted curve shows the full 
contribution up to Q\Z%. The solid curve shows the results of the numerical integration. 

• Figure 2. Same as Figure 1 except for fcsT e = lOkeV. The dotted curve shows the contribution up to Q\Z^. The dashed 
curve shows the contribution up to 9° £ Z$. The dash-dotted curve shows the full contribution up to 9 b e Z 6 . The solid curve shows 
the results of the numerical integration. 

• Figure 3. The ratio T as a function of fcsTe for a fixed value of X = 5. The dotted curve shows the contribution up to 
0fZ4. The dashed curve shows the contribution up to 8lZ$. The dash-dotted curve shows the full contribution up to 8^Zq. 
The solid curve shows the results of the numerical integration. 

This paper has been produced using the Royal Astronomical Society/Blackwell Science BTgX style file. 
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